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What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum 
system? Dodd et al. [Phys. Rev. A 65, 04030f (R) (2002)] provided a partial solution to this 
problem in the form of an efficient algorithm to simulate any desired two-body Hamiltonian evolution 
using any fixed two-body entangling JV-qubit Hamiltonian, and local unitaries. We extend this result 
to the case where the component systems have D dimensions. As a consequence we explain how 
universal quantum computation can be performed with any fixed two-body entangling N-qudit 
Hamiltonian, and local unitaries. 
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I. INTRODUCTION 

A fundamental problem of physics is to determine if 
there exist physical systems that are universal, in the 
sense that they can be used to efficiently simulate any 
other system. A candidate for such a universal system 
was proposed in Deutsch's 1985 paper fjj in the form 
of a universal quantum computer pi G\. The purpose 
of this paper is to investigate what physical systems are 
universal for quantum computation. 

The standard model of a quantum computer consists 
of N qubits, prepared in the state |0)® , that can be ma- 
nipulated by a sequence of one- and two-qubit operations, 
and are subsequently measured in the computational ba- 
sis. There are many possible physical implementations of 
this model, and in general it is an interesting problem to 
determine what critical feature or features of a physical 
system enable universal quantum computation. 

In earlier work by Dodd, Nielsen, Bremner and 
Thew Q it was shown that entanglement is a crucial 
physical ingredient for quantum computation. In partic- 
ular, j| showed that the ability to do local unitary oper- 
ations together with any fixed iV-qubit two-body entan- 
gling Hamiltonian may be used to do universal quantum 
computation on those N qubits. 

In this paper we generalize this result to Hamiltonians 
defined on audits, that is, D-dimcnsional quantum sys- 
tems spanned by the states |0), . . . , \D — 1). This is of 
intrinsic interest, and is also noteworthy because of the 
much richer structure revealed in the general proof than 
in the D = 2 case studied in 01 . 



To state our main result more precisely, we expand an 
arbitrary Hamiltonian on TV qudits as 
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where j%, k±, . . . ,jjsr, fejv each run from through D — 1, 
a jiki...j N k N & re complex numbers, and the operators X 
and Z are /^-dimensional generalizations of the familiar 
Pauli operators, to be defined more precisely later. In 
our work we restrict attention to the case of Hamiltoni- 
ans that only include two-body coupling terms, and do 
not allow three- or more-body coupling terms. A two- 
body coupling between a pair of qudits p and q is a term 
in Eq. (P of the form X^Z k ®X l Z m , where neither X^Z k 
nor X l Z m is equal to the identity, so that the term acts 
nontrivially on qudits p and q, and acts as the identity 
on all other qudits. In order to generate arbitrary en- 
tanglement in the system it is necessary that each qudit 
pair (s, t) must be connected in the sense that there are 
coupling terms in Eq. (Q) for each adjacent pair in some 
sequence (s, . . . , t) of qudits in the system. More explic- 
itly, to any two-body Hamiltonian one can associate a 
graph whose vertices correspond to the qudits in the sys- 
tem, and whose edges connect vertices representing qu- 
dits that are coupled by the Hamiltonian. A Hamiltonian 
is said to be a two-body N-qudit entangling Hamiltonian 
if the graph is connected, that is, there is a path between 
any pair of vertices. Our main result is as follows: 

Let H be a given two-body entangling Hamil- 
tonian on N qudits, and let K be a desired 
two-body Hamiltonian on N qudits. Then we 
have an efficient algorithm to simulate evolu- 
tion due to K using only (a) the ability to 
evolve according to H, and (b) the ability to 
perform local unitaries (that is, single-qudit 
unitaries) on the individual qudits. 
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The algorithm we explain below for performing this 
simulation is only "efficient" in the sense of computer 
science. That is, it requires resources polynomial in the 
number of qudits N in the system. Our simulation tech- 
nique is quite involved, and probably too complicated 
to be experimentally practical. However, the point of 
principle demonstrated by our simulation technique is of 
great importance, namely, that all two-body A^-qudit en- 
tangling Hamiltonians are qualitatively equivalent, given 
the ability to perform local unitary operations. Thus, in 
some sense, the ability to entangle can be regarded as 
a fundamental physical resource — a type of "dynamic 
entanglement" |37j — that can be utilized to perform in- 
teresting processes. We explore this idea and make some 
concrete suggestions for its development in the conclud- 
ing section, Section Furthermore, our work may moti- 
vate future research on more practically viable methods 
for doing universal simulation. 

Antecedents to our work may be identified in many 
different parts of the scientific literature. We now enu- 
merate the different fields in which antecedents may be 
found, before giving a detailed account of the prior work, 
and how it relates to our own. The basic techniques 
we use are generalizations of standard techniques from 
nuclear magnetic resonance (NMR), especially the tech- 
niques known as refocusing and decoupling. The main 
motivation for our work is inspired by research into uni- 
versal gates for quantum computation. More recently 
there has been substantial interest within the quantum 
information science literature in the problem of determin- 
ing when one set of Hamiltonians can be used to simulate 
another. This interest has arisen largely independently 
of work in the quantum control literature, where closely 
related issues are being addressed, albeit using different 
techniques and language. 

The main antecedents of our methods are standard 
NMR techniques for decoupling and refocusing [|[ ||] 
that have been developed and refined over the past half 
century. These techniques have mostly been applied to 
manipulate specific Hamiltonians, rather than general 
classes of Hamiltonians. Ideas from NMR have been ap- 
plied in the quantum computing context by Jones and 
Knill and by Leung, Chuang, Yamaguchi and Ya- 
mamoto [|j , who considered the problem of efficiently im- 
plementing logic gates using a restricted class of Hamil- 
tonians that arises naturally in NMR. 

One of the main motivations for our work is the desire 
to understand what resources are universal for quantum 
computation. Much prior work has been done on this 
subject, and many universal sets of gates for quantum 
computation are known. See, for example [§, The 
work most closely related to ours is independent work of 
J.-L. Brylinski and R. Brylinski jl0|], who used the repre- 
sentation theory of compact Lie groups and real algebraic 
geometry to study the problem of which two-qudit gates 
are universal for computation, given the ability to do 
single-qudit gates. In particular, they defined the class 
of imprimitive unitary gates on two qudits to be the gates 



that are not of the form V ® W or (V ® W)S, where V 
and W are single-qudit unitaries, and S is the swap op- 
eration. They showed that any imprimitive gate is uni- 
versal for quantum computation, given the ability to also 
do arbitrary local unitary operations. Their results thus 
imply ours for the case when N ~ 2. Our results differ 
from theirs in several ways. First, even in the case when 
N = 2, the techniques used in our proof are radically 
different. Our techniques are much more elementary, re- 
lying only on basic linear algebra, a simple result from the 
theory of majorization, and some very elementary num- 
ber theory. Thus, our methods give different insights into 
the problem of universality than those in |u|. Second, 
we consider the case where N > 2, which is potentially of 
great interest for applications to quantum computation 
and quantum control. 

Also related to our results is the work on univer- 
sal gates by Dcutsch, Barenco and Ekert jllj and by 
Lloyd |l2] ], where it was shown that almost any two- 
qubit gate is universal for quantum computation. Lloyd 
sketched a generalization of these results to the case of 
qudits, and this sketch has recently been made rigor- 
ous by Weaver |l3). This work differs from ours in that 
it focuses on unitary gates rather than continuous-time 
Hamiltonian evolution, and does not result in an explicit 
characterization of which sets of unitary gates are uni- 
versal. Our work explicitly determines which two-body 
Hamiltonians, together with local unitary operations, are 
universal. Furthermore, in p| |f2] , 13 it was assumed 
that gates could be independently applied to any pair of 
qudits in the computer, and thus required the ability to 
turn on and turn off interactions between arbitrary pairs 
of qudits. By contrast, we assume only a fixed entangling 
operation. 

Interest in universal quantum gates has recently mo- 
tivated interest in the quantum information science lit- 
erature in the problem of simulating one Hamiltonian 
with another. Independently of Dodd, Nielsen, Bremner 
and Thew M , the problem of Hamiltonian simulation for 
qubits was considered by Diir, Vidal, Cirac, Linden, and 
Popescu in Jl4[ , where it was shown that all two-qubit en- 
tangling Hamiltonians are qualitatively equivalent, in the 
sense that one can be used to simulate the other, given 
the ability to do local unitaries. Wocjan, Janzing and 
Beth Jl5[ considered a specific Hamiltonian acting on a 
system containing A^ spin- 1/2 particles, and considered 
the overhead incurred when using this Hamiltonian to 
simulate other Hamiltonians. Bennett, Cirac, Leifer, Le- 
ung, Linden, Popescu and Vidal | fl6| have considered the 
problem of optimal simulation of one two-qudit Hamil- 
tonian by another, using general local operations, pos- 
sibly including ancillas and measurements. Thus they 
considered a different model than ours, which only in- 
volves local unitary operations, and, in particular, does 
not require the ability to perform interactions with local 
ancilla. Bennett et al showed that in the two-qubit case 
the two models are, in fact, equivalent. We also note that 
the results in [fill are restricted to the case N = 2. Vidal 
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and Cirac |17| extended the results of Bennett et al by 
explicitly obtaining the optimal simulation of one two- 
qubit Hamiltonian by another in the case where classical 
communication between parties is allowed, in addition to 
the ability to do local operations, including the use of 
ancillas and measurement. They also showed that in the 
case of two qudits the model where local unitary oper- 
ations are allowed is distinct from the case where local 
unitary operations and ancilla are allowed, in the sense 
that the latter may be more efficient than the former. We 
note that Leung and Vidal jl8| have independently ob- 
tained results on problems related to those we consider. 
It is worth noting that while related problems are being 
addressed in this long list of papers, the methods used 
are quite varied, and the different methods may provide 
different insights into quantum dynamics. 

Independently of the quantum information science lit- 
erature there has been much interest in Hamiltonian sim- 
ulation in the field of quantum control. A recent overview 
of work in quantum control may be found in |l9[ ] . Of par- 
ticular interest in this context is a general set of neces- 
sary and sufficient conditions for determining whether a 
given set of Hamiltonians can be used to simulate an arbi- 
trary Hamiltonian (see, for example, Schirmer, Solomon 
and Leahy |poft ). These conditions can be applied to de- 
termine whether, in any specific instance, a collection of 
Hamiltonians can be used to simulate an arbitrary Hamil- 
tonian, however they do not directly speak to the ques- 
tion of what class of interactions is universal for quantum 
computation, given the ability to perform local unitaries. 

Finally, we note that the techniques used in this paper 
are closely related to the interesting problem of using a 
Hamiltonian H to simulate time-reversed evolution due 
to the Hamiltonian — H. Results on this problem have 
been obtained by Janzing, Wocjan and Beth and by 
Leung ||. 

The structure of our paper is as follows. Section || in- 
troduces background techniques needed in the main body 
of the paper, including results on the Pauli group and 
Section III explains the N — 2 case of the 
how any two-qudit entangling 



majorization. 
general problem, that is, 
Hamiltonian can be used to simulate any other two-qudit 
Hamiltonian, provided local unitaries are allowed. In Sec- 
tion |y] we explain how this result can be applied to the 
general problem of quantum computation on N qudits, 
and prove our central result. Finally, Section ^ concludes 
the paper with a summary and discussion of our results, 
and a discussion of open problems. 

As the main body of the paper involves a quite exten- 
sive construction, some readers may not wish to wade 
through all the details. We have structured the paper so 
that such a reader may follow the summaries provided at 
the beginning of Sections || and HI , all of Section IV on 
universal quantum computation, and all of the discussion 
m Section [y|. 



II. BACKGROUND 

We now review the background needed to appreciate 
the main body of the paper. At a first read it may be 
useful to skip over the proofs, and pause mainly to appre- 
ciate the nomenclature and basic results. Readers who 
wish to skip the entire section should note the main re- 
sult: given the ability to evolve according to a Hamilto- 
nian J and perform unitary operations Uk it is possible 
to simulate evolution according to a Hamiltonian of the 
form J2k a kUkJUl, where the au are real numbers. This 
composition law for Hamiltonians is the basis for all our 
later simulation results. Note also that throughout the 
paper we use =d to indicate equality modulo D. So, for 
example, 7 =4 3, since 7 is equal to 3, modulo 4. 

The struc ture of the section is as follows. In Sub- 
section [I A we summarize the relations s atisfi ed by op- 
erators in the Pauli group. Subsection [IB describes 



the composition laws used later in the paper to build 
up a library of Hamiltonians we can efficiently simu- 
late given the primitive H amilt onians initially at our dis- 
posal. Finally, Subsection |II C| reviews the basic elements 
of the theory of majorization, including a corollary of 
Uhlmann's theorem crucial to our later analysis. 



A. The Pauli group 

The D-dimcnsional Pauli group consists of all D- 
dimcnsional operators of the form lu XiZ k , where 
j, k, I = 0, . . . , D — 1, ui = exp(27ri / D) , 



X\z) = \z@l)- Z\z)=oj z \z), 



(2) 



and © denotes addition modulo D. The properties of 
the Pauli group were investigated in detail by Gottes- 
man f23|j , and the reader is referred to that paper for 
additional information. 

It is worth noting a few simple properties of the Pauli 
matrices. First, X D = Z D = I. Thus, when writing the 
Pauli matrices we can freely interchange expressions like 
Z D ~ l and Z -1 , and expressions like Z* and Z^ 1 . In a 
similar vein, note that (X j Z k )^ = Z~ k X~K Through 
most of the paper we use notation like Z~ k X~i in pref- 
erence to (x^z k y. 

The basic commutation relations for the Pauli group 
may be written 



(X j Z k ) {X s Z l ) = uj ks ~' Jt (X s z l ) (x j z k ) 



(3) 



We will have very frequent occasion to use these commu- 
tation relations. In particular, note that A- 7 Z k commutes 
with X s Z l if and only if ks =d jt. 

Gottesman Q studied the properties of the Pauli 
group under conjugation by .D-dimensional unitary op- 
erators. In particular, he was interested in normalizer 
operations, that is, unitary operations U such that un- 
der conjugation by U the Pauli group is taken to itself. 
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In Appendix ^| we explicitly describe unitary operators 
performing the following three conjugation operations: 

X -> Z, Z^X- 1 - (4) 
X -> XZ, Z -> Z; (5) 
X->X a , Z -> Z a ~\ when gcd(a, L>) = 1, (6) 



where A — > B means that UAW 



jo 



B for some phase 



factor e l6 . Note that the phase factors are unimportant 
for the proof, and will mostly be ignored in the sequel. 
These equations imply that the following three conjuga- 
tion operations may also be performed: 

X -> Z-\ Z^X; (7) 
X -> X, Z -> XZ; (8) 
X -> X a ~\ Z -> Z a , when gcd(a, D) = 1. (9) 

We now use the normalizer operations to prove what 
we term the Pauli- Euclid- Gottesman (PEG) lemma. 
Aside from its interest as applied in this paper, the PEG 
lemma is also interesting because it enables us to ex- 
plicitly calculate the eigenvalues and eigenvectors of all 
elements of the Pauli group, showing a surprising con- 
nection between the Pauli group and Euclid's algorithm 
(PH, Book 7, Propositions 1 and 2) for finding the great- 
est common divisor. 

Pauli-Euclid-Gottesman Lemma: For 

any dimension D and for integers j and k 
such that 1 < j, k < D — 1, there exists a uni- 
tary operator U such that X j Z k -> Z^ cd ^' k ) 
under conjugation by U . 

Note, incidentally, that the PEG lemma implies that 
the eigenvalues of X^ Z k are equal to the eigenvalues of 
Z scd ^' k \ up to a global phase which may be calculated 
from the proof, below. The eigenvalues of Z^^'^ are 
easily calculable, since Z is already diagonal. The eigen- 
vectors of X^Z may also be extracted from the proof of 
the PEG lemma, below, where we explain how to con- 
struct the conjugating operation, U. 

Froor: From Eqns. (|)-(|) we see that it is possible 
to perform the following operations under conjugation 



X 3 Z k 

x j z k 



X 3 z k+a i 



X 3 



-\-ak iyk 



(10) 

(11) 



where a is any integer. The basic idea of the proof is to 
use these two operations and Euclid's algorithm on the 
paired exponents of X and Z. We will give an example 
of how this is done, with the general proof following sim- 
ilar lines. Consider the operator X 104 Z 80 . Recall how 
Euclid's algorithm is used to find the greatest common 
divisor of 104 and 80. We write 104 = 1 x 80 + 24, so 
gcd(104, 80) = gcd(80, 24). Next, we write 80 = 3x24+8, 
so gcd(104, 80) = gcd(24, 8). Finally, we write 24 = 3 x 8, 
so gcd(104, 80) = 8. These steps are easily mimicked with 



the Pauli operators using Eqns. (|10|) and (|1 1|) . We have 

^-104^80 _^ ^104-1x80^80 _ j^-24^80 ^2) 
X 2i z m ^ X 24 Z 80-3X24 = X 24 Z 8 ^ 

X 2A Z 8 -» X 2i - 3x8 Z 8 = Z 8 . (14) 

The general proof proceeds analogously, using Euclid's 
algorithm. QED 

A key tool in our analysis is the operator expansion. 
We will explain in detail how this expansion works for 
the case of two qudits. Any operator J on two qudits 
may be expanded in the form 



J =z2 r 3k i m X^Z k ®X l Z r ' 

jklm 



(15) 



where the sum is over the range j, k, I, m = 0, . . . , D — 1, 
and the coefficients rjkim may be calculated using the 
expression 



fjklm — 



tr[(Z- k X-3 ® Z- m X- 1 ) J] 



L> 2 



(16) 



In general it is useful to introduce the convention that 
the indices in sums always range over 0, . . . , D — 1, unless 
otherwise noted. 

Eq. ( |l5| ) applies for any operator, however, Hermitian 
operators satisfy additional constraints on the form of 
the coefficients r^im- For example, if a term of the form 
aZ k (g> Z m appears in the operator expansion, then its 
Hermitian conjugate a*Z~ k <X> Z~ m must also appear in 
the operator expansion. In general, terms in the operator 
expansion of a Hermitian operator appear in Hermitian 
conjugate pairs. 

The operator expansion may be used to establish a 
useful identity satisfied by any operator J on a single 
qudit, 



(X j Z k )j(Z- k X-J) =Dtr(J)I. 



(17) 



jk 



The identity Eq. ( |l7| ) is well-known in quantum infor- 
mation science from the properties of the depolarizing 
channel for Z3-dimcnsional systems. The identity may be 
verified by direct calculation, or by substituting an oper- 
ator expansion for J. Eq. ( |T~7| ) may also be extended to 
multiple qudits. For our purposes all that matters is the 
two-qudit case, which reads: 



J2 {x j z k ® x l z m ) j {z- k x- J <g> z- m x- 

jklm 

= D 2 tx{J)I®I. 



(18) 



We conclude this subsection with a brief digression, 
noting that a beautiful alternate proof of Eq. ([IT]) may 
be obtained by applying Schur's lemma from group rep- 
resentation theory f25|| . Let Go denote the Pauli group 
in D dimensions, and note that 



E 

jk 



(x j z k ) j (z- k x- j ) = 1 UJU ^- ( 19 ) 



ueG L 
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The factor 1/D on the right-hand side arises because of 
the phases uo l in front of a general member of the Pauli 
group, uj l X j Z k . The right-hand side of this equation 
commutes with any U G Gd- The result follows from 
Schur's lemma if we can prove that Gd is irreducible. 
Suppose Gd is reducible, so that there exists a non-trivial 
subspace of the qudit state space stable under the op- 
erations in Gd- Let P denote the projector onto that 
subspace. Because the subspace is stable it follows that 
ZPZ~ X = P, and thus Z commutes with P. It follows 
that P can be diagonalized in the basis |0), . . . , \D — 1), 
and thus the stable subspace is spanned by a strict subset 
of |0), . . . , \D — 1). Suppose \z) is in the stable subspace, 
but \z (B 1) is not. But X\z) = \z® 1), so the subspace is 
not stable, which is a contradiction. This completes the 
proof of Eq. (fil). 



(0,0,0,0). Physically, the term -D 2 tr(J)I ® I is an 
unimportant rescaling of the energy and can be neglected. 
The other terms in the expansion for —J are all eas- 
ily simulated using techniques (A) and (B) above. Note 
that the complexity of the simulation scales as 0(D 4 ). 

The above observations (A)-(D) may be summarized 
in a single equation as follows. Given the ability to per- 
form evolution according to the Hamiltonian J and the 
ability to perform unitaries Uk, it is possible to simulate 
evolution according to a Hamiltonian of the form 



(22) 



where the ak can be arbitrary real numbers. 



B. Composition laws for Hamiltonian simulation 

The basic idea employed in the main part of the paper 
is to use our primitive set of operations and a small num- 
ber of composition laws to build up a library of Hamil- 
tonian evolutions we can simulate. We now explain the 
composition laws that we use, adapting from [||. 

(A) Imagine we can evolve according to the Hamilto- 
nian J, and perform unitary operations U and UK Then 
it follows from the identity e _ltc/JC/t = Ue~' ltJ U' { that we 
can exactly simulate evolution according to the Hamilto- 
nian UJlP. 

(B) Imagine we can evolve according to the Hamilto- 
nians J\ and J2. Then we can simulate evolution due 
to J\ + J2 for small times A, due to the approximate 
identity 



-iA(J 1 + J 2 ) 



-iAJi 



-iAJ 2 



(20) 



We treat this identity as though it is exact. This is jus- 
tified, since to simulate a Hamiltonian for a time t it suf- 
fices to perform n separate simulations of a time A = t/n 
each, giving an error of n x 0(A 2 ) = 0(tA). This error 
can thus be made arbitrarily small by making A suffi- 
ciently small. Further remarks on the error analysis are 
made for the qubit case in and these results carry 
over directly to the qudit case. 

(C) Imagine we can evolve according to the Hamilto- 
nian J. Then, by appropriate timing, we can exactly 
simulate evolution according to AJ for any A > 0. 

(D) Imagine we can evolve according to the Hamilto- 
nian J. Then we can evolve according to the Hamilto- 
nian — J. We will explicitly prove this only for the case 
of two-qudit Hamiltonians J, however the proof easily 
generalizes. Note that we can rewrite Eq. dlq), 



j = ( xlzk ® x l z m ) j(z- k x- ] ® z- m x~ 



-Dht(j)i®i, 



(21) 



where we have extracted the — J by taking the sum on 
the right-hand side over all terms except [j,k,l,m) = 



C. Majorization 

The final area of background we shall need is the the- 
ory of majorization, whose basic elements we now review, 
following |2(|. More detailed introductions to majoriza- 
tion may be found in Chapters 2 and 3 of (^7), f2q| , and 
|||. Suppose x = (xi, . . . ,xd) and y = (yi, . . . ,y D ) are 
two D-dimensional real vectors. The relation u x is ma- 
jorized by y" is intended to capture the notion that x is 
more mixed (i.e. disordered) than y. To make the formal 
definition we introduce the notation j to denote the com- 
ponents of a vector rearranged into non-increasing order, 



so x 



to, 



where x\ > x l 2 > 



J 



> x D . We say 



that x is majorized by y, and write x -< y, if 



E 



(23) 



for k — 1, . . . , D — 1, and with the inequality holding with 
equality when k — D. 

The notion of majorization can be extended in a natu- 
ral way to Hermitian operators. We say that the Hermi- 
tian operator A is majorized by the Hermitian operator 
B, and write A -< B, if the spectrum A (A) of eigenvalues 
of A is majorized by the spectrum X(B) of eigenvalues of 
B, where we regard the spectra X(A) and X(B) as vectors. 
So, for example, 



-< 



1 1 

1 ? 

2 2 



(24) 



since the spectra of the two matrices satisfy the majoriza- 
tion criterion, (1/2, 1/2) -< (1,0). 

It is not immediately obvious how this definition of 
operator majorization connects to any natural notion of 
comparative disorder. There is a beautiful result due to 
Uhlmann (30| (see also the reviews J3lj and p2|) that 
provides such a connection. Uhlmann's theorem states 
that A -< B if and only if A = J2 n P n U n BUl, where the 
U n are unitary operators, and the p n form a probability 
distribution. That is, A -< B if and only if A can be 
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obtained from B by mixing together operators unitarily 
equivalent to B. Two important points about the proof 
of Uhlmann's theorem are that the procedure for finding 
the p n and U n is constructive, and, furthermore, there 
are at most D 2 operators U n . 

We now observe that Uhlmann's theorem has a beauti- 
ful corollary when a ppl ied to any two traceless Hcrmitian 
operators A and j?:[p8| 

Theorem: Let A and B be any two trace- 
less Hermitian operators, and assume 6/0. 
Then A -< cB for some positive constant c. 
Uhlmann's theorem then gives an algorithm 
to find a set of at most D 2 unitary operators 
U n , and c„ > 0, such that 



A = J2cnU n BUl 



(25) 



As an example of this theorem in action, consider that 
for any (j, k) ^ (0,0) and (l,m) ^ (0,0) there exist U n 
and c„ such that: 



x?z k + z- k x~i 

Yc n U n (X l Z m + Z- 



l x- 1 ) ut 



(26) 



Using the techniques of the previous subsection, notably 
Eq. (p2|), we see that this equation can be interpreted as 
providing a means of simulating the Hamiltonian X^ Z k + 
Z- k X~i given the Hamiltonian X l Z m + Z- m X- 1 and 
the ability to perform the unitary operations U n . On 
its own, this is not an especially useful simulation result! 
However, similar but more sophisticated variants on this 
idea will be used in our later construction. 

Proof: The case when A = follows by noting that 
-< cB for all c > 0, so we assume A ^ 0. We aim to 
show that X(A) -< c\(B), for some c > 0. Choose 



=i A) (A) 



max 

! 1 " 1 £ ; A;(/i) 



(27) 



Since A and B are traceless and not equal to 0, it follows 
that c > 0. For k = 1, . . . , D — 1 we have 



3=1 



Finally, note that 

D 



r; i -v. ; * \ 

k 

< c£Aj(B). 

3=1 



(28) 
(29) 



£>j(A) = = 

3=1 3=1 

which completes the proof. QED 



(30) 



III. TWO-QUDIT HAMILTONIAN 
SIMULATION 



In this section we study universal simulation with two- 
qudit Hamiltonians, that is, Hamiltonians of the form 



H 



jklrn 



X j Z k ®X l Z r ' 



(31) 



We show that provided this Hamiltonian has a non-zero 
coupling term, that is, a term not of the form / ® (•) or 
(•) (g) /, then H and local unitary operations can be used 
to simulate any other two-qudit Hamiltonian K. 

The basic idea of th e proof is to use the composition 
laws of Subsection [I B to increase the library of Hamilto- 
nians that can be simulated. It will be convenient to use 
the notation Hi, H2, ... to denote the different Hamilto- 
nians that we show how to simulate. The construction is 
rather complicated, for which reason we break it up into 
steps. This separation into steps makes it convenient 
to introduce some global notational conventions. Terms 
like j, k, I, m, n, r, s, t are specific to each step, and some- 
times to individual lines in the proof, often being used 
as dummy variables, with the meaning to be determined 
from context. Terms like a, b, c, d, f carry over from one 
step to another. All of these terms (j, k, . . . and a,b, . . .) 
are consistently integers in the range 0, ... ,13 — 1. 

The general strategy through most of the proof is 
to gradually eliminate more and more terms from the 
Hamiltonian, while keeping particular desired couplings. 
At the end of the proof we are able to simulate a Hamil- 
tonian of an especially simple form, which can then be 
used to build up arbitrarily complicated Hamiltonians. 
We now give an outline of the proof. Note that the 
numbering scheme used in the outline is mirrored in the 
numbering scheme used in the detailed explanation of the 
proof given below in Subsection [II A . 



1. We show that H and local unitaries can be used to 
simulate a Hamiltonian Hi that contains a Z a ® 
Z b coupling term. This term is the focus of most 
of the remaining steps of the proof, as we try to 
eliminate most of the other coupling terms from 
the Hamiltonian. 

2. We show that Hi and local unitaries can be used 
to eliminate terms in Hi not of the form Z' J ® Z k , 
and thus to simulate a Hamiltonian H2 of the form 
J2j k ctjkZi (g) Z k , which still contains the non-zero 
coupling a a bZ a (g) Z b . 

3. We show that H2 and local unitaries can be used 
to simulate a Hamiltonian of the form H3 = 
J2 n Pn(Z c ® Z d ) n , with at least one non-zero cou- 
pling coefficient /3/. 

4. We show that H3 and local unitaries can be used 
to simulate a Hamiltonian of the form Hq = nZ a £g) 
Z b + K*Z~ a ® Z~ b , for any complex n. 
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5. We show that Hq and local unitaries can be used 
to simulate H s = (Z a + Z~ a ) ® (Z b + Z- b ). 

6. Using the corollary to Uhlmann's theorem we show 
that Hg can be used to simulate any Hamiltonian of 
the form Jcg></', where J and J' are arbitrary trace- 
less Hermitian operators. Any two-qudit Hamilto- 
nian can be expressed as a sum of terms of this 
form, together with local interactions, so we con- 
clude that any two-qudit Hamiltonian can be sim- 
ulated using H and local unitaries. 

This construction is complex, and a detailed efficiency 
analysis is not especially enlightening. Nonetheless, from 
the proof below it follows that the total simulation re- 
quires a number of periods of evolution due to H which 
is polynomial in the dimension D. This can be seen by ex- 
amining each step in the construction and verifying that 
they involve only a summation Yk a k(Uk ® Vk)J(Uk <E> 
Vis)' over at most polynomially many terms in D, where 
the Uk and Vk are local unitaries, J is some entangling 
Hamiltonian that we are already able to simulate, and 
the coefficients a.k are also polynomial in D. 

It is worth noting that the proof can be substantially 
simplified if one assumes that D is prime. The reason 
this simplification occurs is because in prime dimensions 
all non-trivial elements of the Pauli group are equivalent 
to one another by unitary conjugation. Thus, given a 
non-zero coupling term in the Hamiltonian it is easy to 
simulate a Hamiltonian in which a non-zero coupling of 
the form Z®Z appears. Given this, steps 1 and 3 can be 
considerably simplified. We descri be in detail how this 
simplification occurs in Subsection III B. 



A. Detailed proof 

Itonian w 
coupling 



1. Simulating a Hamiltonian with a non-zero Z a ® Z b 



By assumption, our Hamiltonian includes a coupling 
term of the form X 3 Z h ® X l Z m with a non-zero coeffi- 
cient. If j ^ or / ^ then the PEG lemma and Eq. (@) 
imply that by performing local unitary conjugation we 
can convert this to a coupling term of the form Z a <& Z b 
with a non-zero coefficient. Let Hi be the Hamiltonian 
that results when this conjugation is performed. It will 
be convenient for our later discussion to fix coprime c and 
d such that c/d = a/b, that is, c/d is a/b in lowest com- 
mon terms. It will also be convenient to define / such 
that a = fc and d = fb. 



2. Simulating a Hamiltonian of the form y]j k ctjkZ-' ® Z k 

We have shown that Hi contains a coupling term of 
the form Z a (& Z b , however it could also contain many 
other coupling terms. We aim to eliminate these other 
terms, while keeping the coupling Z a (£)Z b . In particular, 



we now explain how to eliminate those terms containing 
X or a power of X. Note that given the ability to do Hi 
we can simulate 



H 2 = ^2(Z l ® Z m )Hi(Z-< ® Z~ 



(32) 



To evaluate this sum we first use the commutation rela- 
tions for the generalized Pauli operators and then the ob- 
servation that Y^t wS * ~ D when s —rj 0, and J2t ^ ~ 
otherwise. Using these facts we see that H 2 has the form 

H 2 =J2a jk Z j ®Z k . (33) 
jk 

The term Z a ® Z b in Hi was non-zero, so a a b 0- 



3. Elimination of all terms not of the form (Z c (g) Z ) n 

The next step of the proof is to eliminate all the terms 
in H 2 which are not powers of (Z c £g> Z d ). Note that we 
know there is at least one non-zero term of this form, the 
term Z a (g> Z b = (Z c <g> Z d y . The key to this is a simple 
number-theoretic lemma: 

Lemma: Suppose gcd(Z,m) = 1. Then 
jm =d kl if and only if there exists n such 
that 



j =d nl, and k —d rim. 



(34) 



We will give two proofs of this lemma. The first is a 
constructive proof that only involves elementary number 
theory. The second proof is in some sense more elegant, 
in that it invokes the PEG lemma, and makes use of 
notions of linear algebra. The second proof is given in 
Appendix |b| 

Proof: The reverse implication follows by a simple 
substitution, so we prove only the forward implication. 
Since gcd(Z, m) = 1 there exist integers r and s such that 
rl + sm = 1. Now choose n = jr + ks. Then we have 



nl =£> jrl + ksl 

=d jrl + jsm 

=d j {rl + sm) 

=d j, 



(35) 
(36) 
(37) 
(38) 



as required. A similar calculation shows that nm =d k. 
QED 

Applying our composition laws we see that we can sim- 
ulate 

H 3 = J2{X- d ® X c ) l H 2 {X d ® X- c ) 1 . (39) 



Applying the commutation relations for the Pauli matri- 
ces this simplifies to 



H 3 = a jk 
jk 



(dj-ck)l 



L I 



Z 1 ®Z k . 



(40) 



8 



Note that the sum over I is zero unless dj =n ck. By 
the lemma, this is the case if and only if j =u nc and 
k = D nd for some n, and thus Z j ® Z k — (Z c (g> Z d ) n for 
some n. Thus H% has the form 

H ?> = Y J P n {Z c ®Z d y\ (41) 

n 

where the (3 n are complex numbers. Recall that a = fc 
and b = fd, so (3f oc a a b, and there is at least one non- 
zero coupling term in P 3 . 



4- Simplifying to a sum of at most two terms 

Our next task is to eliminate nearly all the coupling 
terms in H3. First, we set up some notation. Since 
gcd(c, d) — 1 we can choose I and m such that Ic+md = 1. 
It will be convenient to write the coefficients f3 n as a D- 
dimensional vector, that is (3 = ([3q, (3\, . . . , /3d_i), where 
we use the convention that expressions like (x,y,z, . . .) 
denote column vectors. It will also be convenient to use 
the notation e*o, . . . , &d—i f° r the unit vectors in this D- 
dimensional vector space, and to identify e_j with e-D-j- 
So, for example, e\ = (0, 1, 0, 0, ... , 0), and e_2 = e*D-2- 
Note that the constraint that H3 is Hermitian implies 
that /3* = P(3, where Pe n = e_„ for n = 0, . . . , D — 1. 

Next, suppose 7 = (70, . . . , 7.D-1) is a real vector. Us- 
ing our composition laws we can simulate the Hamilto- 
nian 

H 4 =J2lj( x ~ l ® X~ m ) j H 3 {X l ® X rn ) j . (42) 

j 

Our strategy will be to choose 7 in such a way that H4 has 
an especially simple form. Applying the commutation 
relations for the Pauli operators gives 

Hi = J2"f3 0nU {lc+md)3n {Z c ® Z d ) n (43) 

jn 

= J2 S n(Z c ®Z d ) n , (44) 

n 

where in the second step we used the fact that Ic + md = 
1, and we define 

<5 n = /3„^<y%, (45) 
j 

The sum on the right-hand side is most conveniently writ- 
ten in matrix form as -M7, where M is the matrix with 
entries M n j = u) n K Up to a constant M is just the 
matrix representation of the discrete Fourier transform, 
which is easily inverted, so we can choose 7 such that 
M7 = ef + e_/. 

Recall that the 7^ in Eq. ( f42j ) must be real in order 
for a simulation of H4 to be possible. We now use a 
symmetry argument to show that this is the case. Note 



that M7 = ef + (i-f, by definition of 7. Since P, e/ and 
e_y are real, 

M7 - e) + S- f = (P(e> + e_ f ))* (46) 

Next, from M7 = e/ + e_/ and P* = P we obtain 

(P(e> + e_ / ))* = PM*r- (47) 

Combining these results we see that PM*j* = Mj. Ob- 
serving that PM* — M we obtain M7* = M 7, and thus 
7 is real. 

Summarizing, we have obtained the ability to simulate 
a Hamiltonian 

p 4 = /?z a ®^ fc + /rz- a ®z- 6 , (48) 

where [3 = f3f, and the fact that f3-f — (3* follows from 
the fact that H4 is Hermitian. Conjugating by X~ x <S> I 
we also obtain the ability to simulate the Hamiltonian 

P 5 = /3uj a Z a ®Z b + (f3uj a )* Z' a ® Z- b . (49) 

However, note that any complex number k can be formed 
from real linear combinations of (3 and (3u a , so by taking 
appropriate real linear combinations of H4 and P5 we see 
that we can simulate any Hamiltonian of the form 

H 6 = nZ a ®Z b + K*Z~ a ® Z~ b . (50) 



5. Simulation of a tensor product Hamiltonian 

Applying Eq. (||) to the second qudit we see that we 
can simulate any Hamiltonian of the form 

H 7 ± = nZ a ® Z ±b + n*Z- a ® Z Tb , (51) 

and it follows by taking linear combinations that we can 
simulate 

H s = (Z a + Z- a ) ® (Z b + Z- b ). (52) 



6. Simulation of any Hamiltonian 

Note that Z a + Z~ a and Z b + Z~ b are non-zero, 
traceless, Hermitian operators, so by the corollary to 
Uhlmann's theorem we can simulate any Hamiltonian of 
the form J ® J 1 , where J and J' are arbitrary trace- 
less, Hermitian operators. The operator expansion im- 
plies that an arbitrary two-qudit Hamiltonian can be 
formed as a real linear combination of such Hamiltoni- 
ans, together with single-qudit terms of the form J ®I or 
I® J' . Thus, with the ability to perform H and local uni- 
tary operations we can simulate an arbitrary two-qudit 
Hamiltonian. 
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B. The case where D is prime 

The proof just given can be substantially simplified 
in the case where D is prime. We now sketch how the 
simplified proof goes. The reason for the simplification 
is that any non-trivial element X 3 Z k of the Pauli group 
is equivalent under conjugation to Z. To see this, note 
that if j ' and k ^ then, using the PEG lemma 
it is possible to conjugate X 3 Z k to Z, up to a phase 
factor, for some I such that 1 < / < D — 1. Similarly, if 
k = then we can conjugate to some such Z l using (|J), 
while if j = then the term is already in this form. In 
turn this may be conjugated to Z using Eq. (|j|), since 
I is co-prime to D when D is prime. It follows that in 
step 1 of the above proof we can show that it is possible to 
simulate a Hamiltonian Hi that contains a Z®Z coupling 
term. Step 2 proceeds exactly as before, and results in 
a Hamiltonian of the form H2 = J2jk a jk^ 3 <8> Z k , such 
that an 7^ 0. 

Step 3 of the preceding proof is substantially simplified. 
In particular, we note that it is possible to simulate the 
Hamiltonian 

H 3 = Y,( xl ®X- l )H 2 {X- 1 ®X l ) (53) 
1 

= ^2a jk cJ l(k - j) Z j ®Z k (54) 
jki 

= D^a J3 {Z®Zy. (55) 

3 

The remainder of the proof can then be completed as 
before. 



IV. APPLICATIONS TO UNIVERSAL 
QUANTUM COMPUTATION 

We have shown that any two-qudit entangling Hamil- 
tonian, together with local unitary operations, may be 
used to simulate any other two-qudit Hamiltonian. We 
now extend this result to the problem of universal quan- 
tum computation on N qudits. In particular, we show 
that any two-body iV-qudit entangling Hamiltonian, to- 
gether with local unitaries, can be used to perform uni- 
versal quantum computation. 

The basic strategy follows the method presented in j| . 
The idea is to reduce the problem to the two-qudit case 
already solved. To do this, we divide the system into a 
principal system P consisting of two qudits coupled by 
the Hamiltonian H of the entire system, and the remain- 
der of the system, denoted S. Our techniques generalize 
the results in |?], |8| , which are themselves generalizations 
of standard techniques from NMR. The basic idea is to 
turn off all the interactions between P and S, and within 
S, leaving only the interactions present in P. We will 
refer to such a suppression of interactions as decoupling. 
The remaining interactions can then be used, as before, 
to simulate arbitrary dynamics on the two qudits in P. 



Thus it is possible to simulate arbitrary dynamics on any 
two qudits coupled by H. Finally, an arbitrary interac- 
tion between qudits s and t may be effected by perform- 
ing a sequence of SWAP gates between the qudits connect- 
ing s and t (in the sense defined in Section I), applying the 
desired interaction, and swapping back. Note that such a 
sequence of swap operations can be performed using the 
method already described for simulating quantum gates. 
In this way we can effect any two-qudit Hamiltonian be- 
tween any pair of qudits in the system, and thus perform 
universal quantum computation. 

The obvious technique for achieving decoupling is to 
eliminate couplings between P and S, and within S, one 
at a time, using techniques along the lines of those used to 
simulate one two-qudit Hamiltonian with another. Un- 
fortunately, this procedure is not efficient, for reasons we 
now explain. As an example, suppose P consists of two 
qudits, labeled A and B, and S consists of two qudits, 
E and F. Then interactions between qudit E and the 
remainder of the system may be effectively turned off by 
simulating the Hamiltonian 



H' = 



Eu UeHUI 



D 2 



(56) 



where U runs over all Pauli matrices X 3 Z and the E 
indicates that U is being applied to the qudit E. While 
we can in principle turn off all interactions in this way, 
the resulting procedure is not efficient. To see this, notice 
that turning off all the couplings to the qudit E required 
a sum over D 2 terms, each a conjugated form of the 
Hamiltonian H . For an iV-qudit system generalizing this 
procedure in the obvious way would require a sum over 
D N terms. The corresponding simulation would there- 
fore have exponential complexity, which is not efficient. 

Fortunately, much more efficient techniques for decou- 
pling can be devised. In this section we explain two such 
techniques. Subsection IV A explains how the decou- 



pling can be performed for a comple tely a rbitrary two- 
qudit Hamiltonian, while Subsection IV B explains how 



the procedure for decoupling can be substantially simpli- 
fied and made more efficient when the Hamiltonian has 
the localized structure found in most physical systems. 



A. The case of arbitrary two-qudit interactions 

Suppose H is an arbitrary two-qudit entangling Hamil- 
tonian. We now explain how to efficiently eliminate all 
couplings between a principal system P and the remain- 
der of the system S, and to eliminate all couplings inter- 
nal to S, while leaving the couplings within P invariant. 
The method is a straightforward generalization of that 
described for qubits by Dodd et al Q . Let U run over all 
Pauli matrices X 3 Z k . Define Us to be the tensor prod- 
uct of identical operators U acting ditwise on the qudits 
in S. We form the Hamiltonian 



u 



(57) 
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and observe that H' leaves the Hamiltonian on P invari- 
ant, but eliminates all coupling terms between P and S, 
and all single-qudit terms acting within S. 

We now explain a recursive construction to eliminate 
all remaining couplings in S. First, we break the block 
S into two blocks So and S\ of approximately equal size. 
We decouple Sq and Si by forming the Hamiltonian 

u 

Next, we break So into two blocks Soo an d Soi °f approx- 
imately equal size, and break 5*1 into two blocks Sio and 
Sn of approximately equal size. We can decouple Soo 
from 5*10, and Sqi from Six in a single step by forming 
the Hamiltonian 

H "' = J)2 E ( U Soo ® U Sw )H" (U Soo ® C/sJ 1 • (59) 

By repeating this blocking procedure |~log 2 (n — 2)] times 
we can complete the decoupling, leaving a sum over 
O(£)2log 2 JV) = Q( N 2i og2 terms invo i vin g tne conjuga- 
tion of H by local unitary operations. Thus we can de- 
couple P from S, leaving only the interaction on system 
P, using a procedure of complexity 0(N 21 ° 6 ^ D ). This 
interaction on system P can then be used to simulate an 
arbitary two-qudit interaction on P, using the techniques 
described in the previous section. 



B. The case of localized two-qudit interactions 

The method just described assumes a general two- 
qudit Hamiltonian H . Of course, the Hamiltonians oc- 
curring in Nature are usually much more constrained. In 
particular, it is very common for Hamiltonians to have 
some sort of localized structure. In this subsection we ex- 
plain how localized structure can be exploited to obtain 
more efficient decoupling schemes than described above 
for the general case. Note that similar constructions in 
the context of NMR were reported in |7|, || . 

Suppose, for example, that H contains only nearest- 
neighbour interactions on a one-dimensional lattice. This 
is obviously a special case, but is a good illustration of 
the ideas used for more general cases. We number the 
qudits 1, 2, . . . , N, and suppose that P contains qudits 1 
and 2, while S contains qudits 3 through N. The case of 
general P and S follows using similar techniques. We can 
split the decoupling up into three steps. In the first step 
we eliminate all couplings between P and S, which can 
be achieved by eliminating all couplings between qudits 2 
and 3. We call the resulting Hamiltonian H'. The second 
step is to eliminate all single-body terms in S. This can 
be done by simulating the Hamiltonian 

H " = w^2 UsH ' u s- ( 6 °) 

u 



We complete the decoupling by simulating 

H'" = -^^(I<g>/<g>f7<g>/<g>f/<g>I® ...) x 
u 

H'Xl®I®U®I®U®I®...y, (61) 

where the conjugation by U is applied to qubits 3, 5,7,.. ., 
which we can easily see turns off all couplings acting be- 
tween qudits in S. 

Thus, we see that for a nearest- neighbour Hamilto- 
nian on a one-dimensional lattice, the decoupling can be 
performed for constant (with respect to N) cost in the 
simulation, as opposed to the 0(N 2l ° S2 D ) cost incurred 
in the case when general interaction terms appear in the 
Hamiltonian. 

This result can easily be generalized. Suppose S can 
be broken up into a partition Si, ... , S m with the prop- 
erty that qudits in one member of the partition Sj only 
couple to qudits outside Sj. To decouple we do the fol- 
lowing. For each element of the partition Sj turn off all 
couplings between Sj and the remainder of the system 
by simulating the Hamiltonian 

//, y^t . a, t i : . (62) 

u 

where Ho = H . It is easy to see that the Hamiltonian 
H m contains no single-body terms from S, no couplings 
between P and S, and all couplings internal to S have 
been eliminated. The total cost of the simulation scales 
as {D 2 ) m = D 2m . This cost can be reduced even fur- 
ther by using a recursive procedure like that described 
for the general two-qudit case, resulting in a scaling of 
0(D 21og 2 m ). 

Many cases of interest can be described in the 
framework just introduced. For example, consider 
an r-dimensional cubic lattice of qudits, with nearest- 
neighbour interactions. There is a natural partitioning 
of this lattice into 2 r different sublattices, as follows. 
First, fix a site in the lattice, and then consider the cu- 
bic sublattice Si generated by stepping 2 lattice spacings 
in every direction. We generate the partition of sublat- 
tices S\, S2, ■ ■ ■ , SV by translating Si one lattice spacing 
in various directions. (We are ignoring boundary condi- 
tions in this discussion; they are easily accommodated, 
or one can imagine that the lattice has periodic boundary 
conditions) . Now remove the qudits in P from whichever 
elements S p and S q of the partition they happened to fall 
into. Notice that qudits in Sj only ever couple out of Sj, 
since the interactions are nearest-neighbour. Thus the 
procedure described above makes it is possible to decou- 
ple P from S using 0(D 2r ) operations. More generally, 
it is not difficult to use such constructions to efficiently 
decouple P and S for any Hamiltonian containing only 
localized interactions. 
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V. DISCUSSION 

We have shown that, given any two-qudit entangling 
Hamiltonian H and local unitaries we can simulate any 
other two-qudit Hamiltonian. This result was then ap- 
plied to obtain universal gate constructions for quantum 
computation. Our results are of interest because they 
show that such universal simulation is possible, in prin- 
ciple. However, the complexity of our construction limits 
the practicality of potential implementations, and should 
encourage the search for more practical methods. 

There are two aspects to the analysis of efficiency for 
our simulations. The first is how they scale with the 
dimension D of the qudits in the system, and the second 
is how they scale with the number N of qudits present 
in the system. The scaling with N is the critical factor, 
while the scaling with D is not so important, since for 
most physical systems of interest D is a constant. We 
have shown that the scaling for simulation of one two- 
qudit Hamiltonian with another is polynomial in D, and 
the scaling with N behaves as 0(N 21 ° S ^ D ). Thus the 
total scaling is 0(poly(-D)A 21og2 D ), which is polynomial 
in both N and D. 

Our results show that all two-body A-qudit entan- 
gling Hamiltonians are qualitatively equivalent, given the 
ability to perform local unitary operations. Thus, in 
some sense the ability to entangle can be regarded as 
a fundamental physical resource — a type of "dynamic 
entanglement" |39j — that can be utilized to perform in- 
teresting processes. It would be extremely interesting to 
develop a detailed quantitative theory of such dynamic 
entanglement. Following the line of research we have pur- 
sued in this paper, some potential questions one might 
attempt to answer in developing such a theory of dynamic 
entanglement include: 

• What is the optimal procedure for simulating one 
Hamiltonian with another? See |l^, [l(| [l7| for pre- 
liminary results in this direction. 

• Can an entangling Hamiltonian defined onaDxfl' 
system, where D ^ D' , be used to perform univer- 
sal simulation on those systems? Note that this 
question has recently been settled in the affirma- 
tive |33|, using methods rather different than 
that in our paper. 

• Our model assumes that the constituent systems 
are of finite dimensionality D. It would be inter- 
esting to determine whether analogous results hold 
in infinite dimensions. 

• Are universal simulation results possible for non- 
unitary processes? For measurement processes? 
Preliminary results in this direction have been ob- 
tained in @ 

• Can we weaken the condition that arbitrary local 
unitary operations be allowed during the simula- 
tion procedure? It would be interesting, for exam- 



ple, if universal simulation could be performed in 
a system where local unitaries are applied homoge- 
neously across the entire system. 

• Our model assumes that only a single Hamiltonian 
is being applied at any given time, namely, either 
the entangling Hamiltonian H , or a local Hamilto- 
nian on a single qudit. In practice, this is not likely 
to be exactly the case. What effect do imperfec- 
tions have? 

• In the theory of entangled state transformation 
there is a crucial distinction between "single-shot" 
manipulation of entangled states, where just a sin- 
gle copy of the state is available, and manipulations 
that are performed in the asymptotic limit where 
a large number of copies of the state are available. 
The results obtained in the present paper are for 
single-shot Hamiltonian simulation; it would be in- 
teresting to obtain results for the asymptotic case 
as well. 
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APPENDIX A: NORMALIZER OPERATIONS 
FOR THE D-DIMENSIONAL PAULI GROUP 

In this appendix we construct the unitary operations U 
used to perform the conjugation operations Eqs. (|J)-(H), 
which we reproduce here for convenience, 

X -» Z, Z^X- 1 ; (Al) 
X -» XZ, Z -» Z] (A2) 

X -> X a , Z -> Z a ~\ provided gcd(a,Z>) = 1. 

(A3) 

Our constructions are based on those of Gottesman p3| , 
however Gottesman's interest was mainly in the case of 
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prime D greater than 2, and his constructions only apply 
for odd values of D. The following constructions apply 
for D both odd and even. 



The conjugation operation for Eq. (Al) is just the D- 
dimcnsional discrete Fourier transform, defined by 



U\j) 



D-l 

£< 

k=0 



: \k) 



(A4) 



A straightforward calc ulat ion shows that UXU^ = Z and 
UZU^ = Z-\ so Eq. @ holds. 

The definition of the conjugation operation for 
Eq. (A2) depends on whether D is odd or even. When 
D is odd we define 



U\j) 



j0-1)/2 



li>- 



(A5) 



A straightforward calcul atio n shows that UXU^ = XZ 
and UZW = Z, so Eq. (|a|) holds for odd D. When D 
is even we define 



U\j) 



3), 



(A6) 



and then check that UXW = w^IZ and UZW = Z, 
so that Eq. (|a|) also holds for even D. 



Finally, the conjugation operation for Eq. (A3) is de- 
fined by 



U\j) = \aj) 



(A7) 



from which it follows that UXW = X a and UZW = 
Z a , wh ic h co mpletes the constructions needed to verify 
Eqs. 



APPENDIX B: SECOND PROOF OF THE 
NUMBER THEORY LEMMA 



In this appendix we provide an alternate proof of the 



number theory lemma used in Subsection III A of the 
paper. Recall the statement of the lemma: 



Lemma: Suppose gcd(7,m) = 1. Then 
jm =d kl if and only if there exists n such 
that 



j =d n h k = d nm. 



(Bl) 



Proof: By the PEG lemma there exists a normal- 
ize!- operation U such that UX l Z m W = Zs cd C> m ) = Z, 
where the equalities hold up to phase factors. Note 
that X l Z m commutes with X^Z k , since jm —n kl, so 
UX j Z k U^ must commute with Z = UX l Z m Ul It fol- 
lows that UXi Z k U^ — Z n , up to a phase factor, for some 
n, and thus 



x°z k = U^Z n U 
= (U^ZU) n 
= (X l Z m ) n 



(B2) 
(B3) 
(B4) 
(B5) 



where, again, the equalities hold up to unimportant phase 
factors. It follows that j —u In and k =d mn, as 
claimed. QED 
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